Periplasmic chitooligosaccharide-binding protein requires a three-domain organization for substrate translocation

Periplasmic solute-binding proteins (SBPs) specific for chitooligosaccharides, (GlcNAc)n (n = 2, 3, 4, 5 and 6), are involved in the uptake of chitinous nutrients and the negative control of chitin signal transduction in Vibrios. Most translocation processes by SBPs across the inner membrane have been explained thus far by two-domain open/closed mechanism. Here we propose three-domain mechanism of the (GlcNAc)n translocation based on experiments using a recombinant VcCBP, SBP specific for (GlcNAc)n from Vibrio cholerae. X-ray crystal structures of unliganded or (GlcNAc)3-liganded VcCBP solved at 1.2–1.6 Å revealed three distinct domains, the Upper1, Upper2 and Lower domains for this protein. Molecular dynamics simulation indicated that the motions of the three domains are independent and that in the (GlcNAc)3-liganded state the Upper2/Lower interface fluctuated more intensively, compared to the Upper1/Lower interface. The Upper1/Lower interface bound two GlcNAc residues tightly, while the Upper2/Lower interface appeared to loosen and release the bound sugar molecule. The three-domain mechanism proposed here was fully supported by binding data obtained by thermal unfolding experiments and ITC, and may be applicable to other translocation systems involving SBPs belonging to the same cluster.

Binding mode of (GlcNAc) 3   Figure 3A shows a 2F 0 -Fc map of the bound (GlcNAc) 3 .As in the case of VhCBP 29 , electron density of (GlcNAc) 3 was identified in the (GlcNAc) 3 -liganded VcCBP.Since Kitaoku et al. 30 designated these subsites as Site1, Site2 and Site3, the same subsite nomenclature was used in this study.In the (GlcNAc) 3 -liganded VcCBP, the electron density of the Site3 GlcNAc was lower than those of the two GlcNAc residues at Site1 and Site2, and the occupancy of Site3 GlcNAc was set at 0.5.Two conformers of the Phe410 side chain were observed near the Site3 GlcNAc, and the individual occupancies of the side chain were also set at 0.5.It is most likely that the mobility of the Site3 GlcNAc is much higher than that of the other GlcNAc residues and that the Phe410 side chain flips alternatively back and forth in the liganded structure.Both GlcNAc units at Site1 and Site2 appeared to be hydrogen-bonded or in hydrophobic contact with Asp9, Asn203, Ser220, Phe221, Trp362, Asp364, Asn408, Phe410, Arg435, Phe436 and Trp512, of which all except Asp9 are conserved between VhCBP and VcCBP (highlighted on a black background in Figs. 1 and 3B), while the GlcNAc unit at Site3 makes only a couple of hydrogen bonds with the main chain carbonyl of Ala513 and the guanidyl nitrogen of Arg27, which are localized to the interface between the Upper1 and Upper2 domains (Fig. 3B, Supplementary Table S1).
B-factor values visualized in the main chain structure of (GlcNAc) 3 -liganded VcCBP are shown in Fig. 4A, and were higher in the Upper2/Lower interface as well as the Site3-GlcNAc contact surface of Upper1 domain.Notably, the B-factors were significantly higher at the loop region immediately following the secondary structure η6 and the loop between α14 and α15 (red broken circles, Fig. 4A).These two loops are located at the Upper2-Lower interface.Furthermore, B-factor values of the individual GlcNAc residues were 9.86, 12.67 and 16.28 Å 2 , respectively (Fig. 4D), indicating the higher mobility of the Site3 GlcNAc.These B-factor data were fully www.nature.com/scientificreports/consistent with the electron density data described above.The cooperative motions of the Site3 GlcNAc and the Phe410 side chain may be significant from the mechanistic viewpoint.
The unliganded form exhibited the largest motion (RMSD, 0.8 nm; the top panel of Fig. 5A) of the entire protein molecule, while the motions of the individual domains remained at lower levels (RMSD, 0.1-0.3nm; Fig. 5A).This clearly indicated a domain motion, in which the individual domains did not change their own conformations but changed their relative arrangement.On the other hand, the motions of the entire protein molecules as well as individual domains were very low in the (GlcNAc) 2 -liganded and (GlcNAc) 3 -liganded VcCBP, where RMSDs were only 0.1-0.2nm (Fig. 5B,C).Figure 6 shows the distances between the α-carbons of the two numbered amino acid residues (amino acids labeled in red, Fig. 3B).The distance 203-361, which reflects the separation between Upper1 and Lower domains (red broken line between Asn203 and Gly361 in Fig. 3B; the top panels of Fig. 6A,B), did not significantly differ between (GlcNAc) 2 -liganded and (GlcNAc) 3 -liganded VcCBP.For the distance 24-391 (red broken line between Thr24 and Ala391 in Fig. 3B; second panels from the top of Fig. 6A,B), which also reflects the separation between Upper1 and Lower domains at another site, no significant difference was observed between (GlcNAc) 2liganded and (GlcNAc) 3 -liganded VcCBP.Furthermore, the distance 135-517 (red broken line between Ser135 and Glu517 in Fig. 3B; the bottom panels of Fig. 6A,B) reflecting the separation between the Upper1 and Upper2 domains, did not exhibit any differences between the two liganded VcCBP.However, the distances 101-432 (Gln101-Gly432) and 24-101 (Thr24-Gln101) reflecting the separations of Upper2-Lower and Upper1-Upper2, respectively, fluctuated more in (GlcNAc) 3 -liganded VcCBP than in (GlcNAc) 2 -liganded (the third and fourth Disallowed region (%) 0 0 panels from the top of Fig. 6A,B).The movements are more intensive in the (GlcNAc) 3 -liganded VcCBP than in the (GlcNAc) 2 -liganded in this region.Elongation of the chain of the bound oligosaccharide from (GlcNAc) 2 to (GlcNAc) 3 was found to enhance the molecular motion around amino acid residue Gln101 located in the Upper2 domain.The larger conformational change in (GlcNAc) 3 -liganded VcCBP was also confirmed from Fig. 7, which shows the 2-dimensional projection of the trajectories of principal components, eigenvector1 and eigenvector2.

Thermal unfolding experiments with VcCBP in the presence of (GlcNAc) n
A typical dataset of the thermal unfolding transitions of VcCBP in the absence or presence of (GlcNAc) n was shown in Supplementary Fig. S1.The unfolding transition was highly cooperative, with the fraction unfolded increasing sharply at the transition temperature (T m ).The individual unfolding experiments were conducted twice, and the averaged values of T m were calculated as listed in Table 2. Thermal unfolding of unliganded VcCBP took place at T m of 43.2 °C.When (GlcNAc) 2 , (GlcNAc) 3 or (GlcNAc) 4 was added to VcCBP, T m was elevated from 43.2 to 54.4 °C (ΔT m = 11.2 °C), 54.8 °C (ΔT m = 11.6 °C) and 53.8 °C (ΔT m = 10.6 °C), respectively.As in the case of VhCBP 29 , these three oligosaccharides bound strongly to VcCBP.However, ΔT m values induced by the addition of (GlcNAc) 5 and (GlcNAc) 6 were moderate, 6.3 °C for both, indicating that these two oligosaccharides bind to VcCBP with significantly lower affinities 34 .ΔT m for GlcNAc was only 2.8 °C, suggesting that the monomer bound to VcCBP with much lower affinity.
Thermodynamic parameters for binding of (GlcNAc) n (n = 2, 3, and 4) ITC experiments were conducted at 25 °C by titrating a 1.0 mM solution of GlcNAc or (GlcNAc) n (n = 2, 3, 4, 5 or 6) into VcCBP solution (50 µM).As shown in Fig. 9A-F, titrations with GlcNAc, (GlcNAc) 5 , or (GlcNAc) 6 did not result in any heat release/absorption, whereas (GlcNAc) 2 , (GlcNAc) 3 and (GlcNAc) 4 exhibited clear heat release or absorption.This was consistent with the experimental data obtained from thermal unfolding experiments (Supplementary Fig. S1 and Table 2).It is noteworthy that titration with (GlcNAc) 2 resulted in heat release but (GlcNAc) 3 and (GlcNAc) 4 resulted in absorption.The thermodynamic mechanism of the interaction with (GlcNAc) 2 appears to be basically different from those with (GlcNAc) 3 and (GlcNAc) 4 .Thermodynamic parameters were obtained based on the experimental data as listed in Table 3.For (GlcNAc) 2 , the enthalpic contribution (ΔH°) and the association constant (K assoc ) were obtained based on the theoretical fit and were − 2.40 kcal mol −1 and 4.08 × 10 6 M −1 , respectively.Thus the binding free energy (ΔG°) was − 8.85 kcal mol −1 .The entropic contribution (− TΔS°) was calculated to be − 6.45 kcal mol −1 .The interaction was entropy-driven with a smaller enthalpy gain.When (GlcNAc) 3 or (GlcNAc) 4 was used instead of (GlcNAc) 2 , the entropy gain was even greater but this was compensated by the loss of enthalpy, resulting in lower binding affinities (ΔG°) of − 7.54 and − 8.36 kcal mol −1 , respectively.3).The results are also listed in Table 4.Both favorable ΔS°s olv and unfavorable ΔS°c onf contributed to the interaction with (GlcNAc) n .The favorable ΔS°s olv can be explained by the exclusion of bound water molecules upon ligand binding, while the unfavorable ΔS°c onf by the reduction in the fluctuation of the VcCBP structure as seen in Fig. 5, which shows the fluctuations of RMSDs in unliganded and liganded states.

Binding experiments using VcCBP_R27A
To define the contribution of the Arg27 side chain to (GlcNAc) 3 binding, we conducted ITC analysis of the interaction with (GlcNAc) 3 using VcCBP_R27A.As in the case of VcCBP, heat absorption was observed up to 20 °C, and converted to heat release at 25-35 °C (Supplementary Fig. S4).The thermodynamic parameters were obtained for the individual temperatures, and the values were plotted as shown in Supplementary Fig. S3D.At all temperatures tested, the binding affinities (ΔG°) were found to be reduced by about 1 kcal/mol in VcCBP-R27A (Supplementary Fig. S3D, Table 3).Thus Arg27 participates in the binding of (GlcNAc) 3 but its contribution is only moderate.The temperature dependence of ΔH° of VcCBP_R27A indicated that both the solvation entropy gain and conformational entropy loss were reduced (ΔS o solv , 110.9 → 93.5 cal K −1 mol −1 ; ΔS o conf , − 70.3 → − 59.9 cal K −1 mol −1 ), and the reduction in the total entropic contribution (− TΔS o , − 9.56 → − 7.30 kcal mol −1 ) was compensated by the reduction in enthalpy loss (ΔH o , 2.02 → 0.92 kcal mol −1 ) (Tables 3 and 4).

Discussion
(GlcNAc) 2 binds primarily to the Upper1/Lower interface (Site1/Site2) As listed in Table 3, the interaction of (GlcNAc) 2 with VcCBP was entropy-driven (− TΔS° = − 6.45 kcal mol −1 ) with a smaller enthalpy gain (ΔH° = − 2.40 kcal•mol −1 ).This is consistent with the data reported for VhCBP (ΔH° = − 3.9 kcal•mol −1 and − TΔS° = − 6.4 kcal mol −1 ) 28 .Since the amino acid residues involved in (GlcNAc) 2 binding in VhCBP are almost conserved in VcCBP as shown in Fig. 1, the mechanism of (GlcNAc) 2 binding in VcCBP is similar to that in VhCBP from the viewpoints of structure and thermodynamics.As shown in Fig. 3B and Supplementary Table S1, only the Upper1/Lower interface appears to be involved in the binding of two GlcNAc residues from the non-reducing end.Since VcCBP has been regarded as being specific for (GlcNAc) 2 translocation to the ABC-transporter 18 , (GlcNAc) 2 binds primarily to the interface between the Upper1 and Lower domains (Site1/Site2) before translocation.

Enthalpic Site1/Site2 and entropic Site3/Site4
From the thermal unfolding and ITC experiments, we found that the binding affinities of (GlcNAc) n (n = 2, 3, and 4) were comparable with each other (Supplementary Fig. S1 and Fig. 9; Tables 2 and 3), but that favorable contributions of entropic term (− TΔS°) to the binding affinities were enhanced in the order (GlcNAc) 2 < (GlcN Ac) 3 < (GlcNAc) 4 (Table 3).As the entropic contributions were enhanced, compensations were clearly found in the enthalpic terms.As seen from the crystal structures of (GlcNAc) 2 -liganded (1ZU0) and (GlcNAc) 3 -liganded VcCBP (Fig. 3B), reducing-end GlcNAc of bound (GlcNAc) 3 exists beyond the Upper1/Lower interface and that of bound (GlcNAc) 4 may be in contact with the Upper2/Lower interface (Site3 and Site 4).From the temperature dependence of ΔH°, we found that favorable solvation entropy (ΔS°s olv ) predominated over conformational entropy changes (ΔS°c onf ) as listed in Table 4, suggesting that Site3 and Site4 bind a relatively large number www.nature.com/scientificreports/ of water molecules in the open form, and that the bound water molecules may be excluded upon (GlcNAc) 3 or (GlcNAc) 4 binding.We calculated the solvent-accessible surface areas (ASA) for unliganded, (GlcNAc) 2liganded, and (GlcNAc) 3 -liganded VcCBPs based on their crystal structures 36 , and the data were presented in Supplementary Table 2.We found that the apolar solvent accessible surface area (ASA apolar ) was reduced by 3.3 or 4.2% upon binding of (GlcNAc) 2 or (GlcNAc) 3 , respectively.This indicated that (GlcNAc) 3 binding to VcCBP excludes more water molecules from the apolar surface than (GlcNAc) 2 binding.This is consistent with the positive values of ΔS solv , which is greater in (GlcNAc) 3 than in (GlcNAc) 2 (Table 4).The thermodynamic data

Upper1 and Upper2 domains play different roles in (GlcNAc) 2 translocation
In VcCBP, we found three structural domains, the Upper1, Upper2 and Lower domains (Fig. 2), which correspond to domains I, II and III of cluster C SBPs reported by Chandravanshi et al. 19 .The three-domain organization appears to be significant from a functional viewpoint.One of the reasons for this significance derives from the clear distinction found in the binding thermodynamics between Site1/Site2 and Site3/Site4; although both enthalpy and entropy contributions were involved in the interaction with the former site, the entropic contribution predominates in the latter site (Tables 3 and 4).The second reason is the larger RMSD of the Upper2 domain upon superimposition of (GlcNAc) 2 -liganded and (GlcNAc) 3 -liganded VcCBP structures (Fig. 2C).The larger RMSD is likely to produce the space for accommodating the additional GlcNAc unit at the interface between Upper1 and Upper2 domains.The third reason is the distinction in the molecular movements between the Upper1/Lower and Upper2/Lower interfaces.The distance Gln101-Gly432 fluctuated more in (GlcNAc) 3liganded than in (GlcNAc) 2 -liganded VcCBP, whereas the distances, Asn203-Gly361 and Thr24-Ala391 were similar (Fig. 6).In (GlcNAc) 3 -liganded VcCBP, the fluctuations were more intensive in the Upper2/Lower interface (Gln101-Gly432) than in the Upper1/Lower interface (Asn203-Gly361/Thr24-Ala391) as shown in Fig. 6B.The larger fluctuation in agreed well with the highest B-factor values in the loop immediately following η6 and in the loop between α14 and α15 (Fig. 4).A "half-open" conformation observed in Trp513-mutated VhCBP 30 may correspond to a snapshot of the largely fluctuated state.Cross-correlation heat maps of (GlcNAc) 2 -liganded and (GlcNAc) 3 -liganded VcCBP (Fig. 8) also showed the independency of the Upper1 and Upper2 domains in their molecular movements; thus, it is most likely that Upper1 and Upper2 domains play different roles in the (GlcNAc) 2 translocation process.

Hypothetical releasing site, Site5/Site6
As shown in Fig. 9, titrations of VcCBP with (GlcNAc) 5 or (GlcNAc) 6 did not result in any heat release/absorption; the binding affinities are too low to obtain the thermodynamic parameters for these oligosaccharides by ITC.Nevertheless, thermal unfolding data (Supplemental Fig. S1 and Table 2) revealed a significant elevation of the transition temperature of thermal unfolding (ΔT m = 6.3 °C), suggesting a significant interaction of (GlcNAc) 5 or (GlcNAc) 6 with VcCBP.In the interaction of (GlcNAc) 6 with VcCBP, the two GlcNAc residues of the non-reducing end interact with the interface between the Upper1/Lower (Site1/Site2) and the neighboring (GlcNAc) 2 unit also interacts with the Upper1/Upper2/Lower interface (Site3/Site4).However, the additional (GlcNAc) 2 unit of the reducing-end side may be repelled from the protein surface (hypothetical Site5/Site6).This situation may Table 3. Thermodynamic parameters obtained for the VcCBP/(GlcNAc) n (n = 2, 3, and 4) interactions by means of ITC.The individual ligands (1.0 mM) were titrated into the VcCBP or VcCBP_R27A solution (50 µM) in 20 mM Tris-HCl buffer, pH 8.0, in the sample cell (0.2028 mL).The experiments were performed at 25 °C using an iTC200 system (Microcal Northampton, MA, USA).Three repeated experiments were conducted to obtain individual sets of the thermodynamic parameters.2, accounting for the lower binding affinity of (GlcNAc) 5 and (GlcNAc) 6 .We propose here that a specific substrate (GlcNAc) 2 primarily binds Site1/Site2 with both enthalpyand entropy-driven interactions, and is subsequently translocated to Site3/Site4, where the binding interaction is looser, leading to release of the sugar molecule from Site5/Site6 to a (GlcNAc) 2 -specific ABC transporter.The loosening of the interaction at Site3/Site4 may be caused by the higher mobility of the Upper2/Lower interface observed in the (GlcNAc) 3 -liganded structure (Figs. 4, 6 and 7).Thus, all translocation processes are conducted by the cooperative action of the three domains, Upper1, Upper2 and Lower.

Structure triggering the (GlcNAc) 2 unit translocation
Kitaoku et al. 30 observed electron density of the Phe411 side chain with a full occupancy (1.0) close to the N-acetyl methyl group of the Site2 GlcNAc in (GlcNAc) 2 -liganded VhCBP.However, in (GlcNAc) 3 -liganded VhCBP, the same side chain was found 3.5 Å away from the Site2 GlcNAc with a full occupancy (1.0).Here, in (GlcNAc) 3 -liganded VcCBP we observed the side chain of the corresponding phenylalanine residue (Phe410) at both positions with individual occupancies of 0.5 (Fig. 3A), indicating a significant flipping of the Phe410 side chain.The density of the Site3 GlcNAc was also observed with an occupancy of 0.5, indicating the higher mobility of the Site3 GlcNAc (Figs. 3A and 4D), which appeared to be coordinated with that of the Phe410 side chain.The coordinated motions of the phenylalanyl side chain and the Site3 GlcNAc suggested that Phe410 may be involved in the translocation of bound (GlcNAc) 2 from Site1/Site2 to Site3/Site4.Perhaps the translocation process is triggered by interaction with the corresponding ABC transporter, which may further translocate (GlcNAc) 2 to the hypothetical release site, Site5/Site6, located in the Upper2/Lower interface.

Conclusion
Taken together, the three domains, Upper1, Upper2, and Lower domains, found in the crystal structure of VcCBP play different roles and function cooperatively in translocation of (GlcNAc) 2 .The mechanism proposed here was fully supported by binding data obtained by thermal unfolding and ITC experiments and may be applicable to other translocation systems involving SBPs belonging to the same cluster.

Construction of expression plasmid for VcCBP and VcCBP_R27A
Synthetic genes encoding VcCBP and VcCBP_R27A, in which Arg27 was mutated to alanine, were obtained from Invitrogen (Carlsabad, CA, USA).The nucleotide sequences of the genes were optimized to increase expression in E. coli without changing the amino acid sequences of these proteins.The expression vectors, pRham-VcCBP and pRham-VcCBP_R27A, were constructed by the Expresso (R) Rhamnose Cloning and Expression System, N-His (Lucigen, UK).

Protein expression and purification
The expression vector, pRham-VcCBP or pRham-VcCBP_R27A was transformed into E. coli C43(DE3).Induction with 0.2% α-L( +)-rhamnose was conducted according to the supplier's instruction.After induction, the culture was incubated at 15 °C for 40 h, then the cells were harvested and disrupted by sonication in 20 mM Tris-HCl buffer, pH 8.0.The sonicated extract was centrifuged at 12,000×g for 15 min at 4 °C.The soluble fraction was dialyzed against 20 mM Tris-HCl buffer pH 8.0, and applied to a Ni-NTA column equilibrated with the same buffer.After washing the column with the Tris buffer, the bound protein fractions were eluted with a linear gradient of 0-0.2 M imidazole.The fractions containing a protein of molecular mass 60 kDa, which corresponds to that of VcCBP, were collected, and ammonium sulfate was added to the protein solution to a final concentration of 1 M.The solution was applied to the TOYOPEARL Butyl 650 M column equilibrated with 20 mM Tris-HCl buffer pH 8.0 containing 1 M ammonium sulfate.The adsorbed fraction was eluted with a linear gradient of 1-0 M ammonium sulfate in 20 mM Tris-HCl buffer pH 8.0.The protein fractions containing VcCBP were pooled and applied to a Q-Sepharose column previously equilibrated with the same buffer.The protein was eluted stepwise with 0.15 M NaCl in the same buffer.The fractions containing VcCBP were pooled and further applied to a Sephacryl S-100 HR gel-filtration column equilibrated with the same buffer containing 0.1 M NaCl.Fractions exhibiting a single protein band on SDS-PAGE 38 (Supplementary Fig. S5) were pooled and stored at 4 °C.

Protein concentration
Protein concentrations were determined by reading absorbance at 280 nm, using an extinction coefficient of VcCBP (110,365 M −1 cm −1 ) calculated from the equation proposed by Pace et al. 39 .

Thermal unfolding experiments
To obtain the thermal unfolding curve of VcCBP, the CD value at 222 nm was monitored using a Jasco J-720 spectropolarimeter (cell length 0.1 cm), while the solution temperature was raised at a rate of 1 °C min −1 using a temperature controller (PTC-423L, Jasco).

Isothermal titration calorimetry (ITC)
Solutions of 50 μM VcCBP or VcCBP_R27A in 20 mM Tris-HCl buffer, pH 8.0, were degassed and loaded into the sample cell (0.2028 mL).The individual ligands (1.0 mM) were dissolved in the same buffer, degassed and loaded into a syringe.Calorimetric titration was performed at 25 °C using an iTC200 system (Microcal Northampton, MA, USA).In the titrations, 2.5 μL of a ligand was injected into the sample cell at 180-s intervals with a stirring speed of 1000 rpm.The dilution heat caused by each titration was measured by titrating ligand to buffer solution without protein under identical conditions.The dilution heat thus obtained was subtracted from the heat change that was observed in the presence of protein.Individual titration experiments were repeated three times to obtain reliable values of thermodynamic parameters.The Origin software installed in the ITC instrument was used to analyze the ITC data.Using the One-set of Sites model, individual datasets obtained from the titration experiments fitted well to the theoretical curves, providing the stoichiometries (n), equilibrium association constants (K a ) and enthalpy changes (ΔH°) of the protein-ligand interactions.The binding free energy change (ΔG°) and entropy change (ΔS°) were calculated from the relationship as follows, The accuracy of the thermodynamic parameters obtained was assessed from the c-values calculated from the equation c = n•K a •[M] t , where [M] t is the total concentration of protein 41 .The ITC measurements were conducted in 20 mM Tris-HCl buffer pH 8.0, at various temperatures from 5 to 35 °C.The ΔH° values obtained for various temperatures were plotted against temperatures, and the slope of a straight line fitted to the experimental points corresponds to the heat capacity change (ΔC p °).As the entropy of solvation is regarded as zero for proteins near 385 K, ΔC p ° was converted to the solvation entropy change (ΔS solv °) at 25 °C (298 K) according to the following relationship, The conformational entropy change (ΔS conf °) was calculated from ΔS° obtained from Eq. (1), the solvation entropy change (ΔS solv °) and the mixing entropy change (ΔS mix °, − 8 cal K −1 mol −1 ), based on the following Eq. 42,

Crystallization and data collection
Crystallization conditions for VcCBP were screened using the sparse-matrix sampling method by sitting drop vapor diffusion at 20 °C.Under optimized crystallization conditions, 1 μL of protein solution (5 mg/ml in water) was mixed with 1 µL of 0.1 M sodium citrate containing 0.2 M ammonium acetate with 30% w/v polyethylene glycol 4000, pH 5.6.Rod-like crystals of VcCBP grew within 3 weeks to a size of up to 0.1 × 0.1 × 0.5 mm 3 .To prepare the crystals of (GlcNAc) 3 -liganded VcCBP, the unliganded crystals were transferred to the crystallization well solution containing 26 mM (GlcNAc) 3 and incubated at 20˚C for 3.5 h.The crystals were successfully grown in the presence of (GlcNAc) 3 .For data collection, the crystals were transferred into the cryoprotectant solution containing 0.01 M zinc sulfate, 0.1 M MES (pH 6.5) and 30% PEG MME550, and then flash-cooled in a nitrogen stream at 95 K.The diffraction data were collected at the beam-line BL-17A of Photon Factory (Ibaraki, Japan), using an EIGER X 16 M (Dectris), at a cryogenic temperature (95 K).The data were integrated and scaled with XDS 43 .The processing statistics are summarized in Table 1.

Structural determination and refinement
The structures of unliganded and (GlcNAc) 3 -liganded VcCBP were solved by the molecular replacement method using the program PHASER 44 , where the structures of unliganded VcCBP (PDB code, 1ZTY) and (GlcNAc) 2liganded VhCBP (PDB code, 5YQW) served as search models, respectively.For unliganded VcCBP, two protein molecules were located in the crystallographic asymmetric unit.The model was improved by several rounds of refinement with PHENIX 45 and COOT programs 46 .Occupancies of the Site3 GlcNAc and two conformers of the Phe410 side chain were set at 0.5.The structure of unliganded VcCBP was refined to an R work /R free of 16.2/19.1% at a resolution of 1.6 Å.The final model contains two protein molecules that include residues 1-532 for each molecule and 773 water molecules.The stereochemistry of the model was verified using MolProbity 47 , showing 96.7%, 3.3% and 0% of protein residues in the favored, allowed and disallowed regions of the Ramachandran plot, respectively.For (GlcNAc) 3 -liganded VcCBP, one protein molecule was located in the crystallographic asymmetric unit.The structure of (GlcNAc) 3 -liganded VcCBP was refined to an R work /R free of 16.9/18.7%at a resolution of 1.22 Å.The final model contained one protein molecule that includes residues 1-532 and 563 water molecules.The stereochemistry verification showed 97.3%, 2.7% and 0% of protein residues in the corresponding regions, respectively.Molecular graphics were illustrated with PyMol software (http:// www.pymol.org/).The refinement statistics are summarized in Table 1.

Figure 3 .
Figure 3. State of (GlcNAc) 3 bound to VcCBP.(A) 2F 0 -F c map of the bound (GlcNAc) 3 and the Phe410 side chain.The occupancies of the Site3 GlcNAc and two conformers of the Phe410 side chain were each assumed to be 0.5.(B) Binding mode.The color system is identical to that in Fig. 2. The possible hydrogen bonds are shown by black broken lines.The red broken lines represent the distances between the α-carbons monitored in the molecular dynamics simulation shown in Fig. 6.

Figure 4 .
Figure 4. Visualization of crystallographic B-factors for the VcCBP structure.(A) and (B) entire VcCBP structure in complex with (GlcNAc) 3 ; (C) Cross-section of the (GlcNAc) 3 -binding site; and (D), close-up view of the (GlcNAc) 3 -binding site.The colors from violet to blue, green, yellow, orange and red indicate B factors from small to large.Broken circles colored white represent B-factor-higher regions, while those colored red represent loop regions with the highest B-factors.

Figure 8 .Table 2 .
Figure 8. Cross correlation heat maps of Cα atoms around their mean positions for the entire simulation period.(A) Unliganded VcCBP (8I5J); (B) (GlcNAc) 2 -liganded VcCBP (1ZU0); and (C) (GlcNAc) 3 -liganded VcCBP (8I5K).The color gradation from red to blue corresponds to extents of correlated motions (from 1 to 0, positive correlations) and anti-correlated motions (from 0 to − 1, negative correlations).The horizontal and vertical axes represent the amino acid residue No. of VcCBP.The thin guided lines in the figures represent the boundaries of the individual domains, Upper1, Upper2, and Lower.The correspondence table on the right shows which domain is located where.
Vol:.(1234567890) Scientific Reports | (2023) 13:20558 | https://doi.org/10.1038/s41598-023-47253-ywww.nature.com/scientificreports/bring about the lower but significant ΔT m in Table 40olar ellipticities) were normalized to obtain unfolded fractions at individual temperatures.To assess the binding ability of GlcNAc and (GlcNAc) n (n = 2, 3, 4, 5, and 6), the unfolding experiments of VcCBP were conducted in the presence or absence of (GlcNAc) n .Individual unfolding experiments were repeated twice under the same conditions.The transition temperature of thermal unfolding (T m ) was elevated when the ligand was added to the VcCBP solution.The elevation of T m (ΔT m ) indicated the binding of the ligand40.The solvent condition was 20 mM Tris-HCl buffer pH 8.0.Final concentrations of the protein and (GlcNAc) n were 8 µM and 8 mM, respectively.